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In this work we study in detail the NeVII 2s3p-2s3s singlet line, which was also the object of a 
recent experiment. The standard perturbative impact theory predictions are tested against a fully 
non-perturbative semiclassical impact calculation, taking into account dipole and quadrupole inter- 
actions. Potentially very significant problems with the standard perturbative theory are encountered 
and discussed and a simple remedy is proposed. 
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The calculation of plasma-broadened line spectra, pro- 
vides a very useful diagnostic tool and additionally is a 
necessary ingredient for large-scale computations in as- 
trophysics and plasma physics. A major cornerstone, re- 
ducing the many-body problem of line-broadening to the 
computation of one-body quantities, is the impact ap- 
proximation [|l],^). For practical calculations, the impact 
theory is usually employed in its perturbative version, 
\Q ' and even then simplified formulas are often used. 

Isolated lines by being relatively simple and usu- 
ally unaffected by ion microfield effects, whether static 
or dynamic || are an excellent testing ground for the 
theories of electron collisional broadening. Such tests re- 
quire reliable experimental profiles and much progress 
has been made in this direction in recent years mainly 
by the Bochum group |i|-|ic|l . These studies have revealed 
significant discrepancies with simplified expressions that 
are often used for the electron collisional broadening 
|@,[EL}jl4). Furthermore, serious discrepancies with close 
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coupling(cc) calculations jL5) were found in ||. Even 
more impressive is a recently obtained factor of 2 discrep- 
ancy between cc calculations jl5| and experiment [ flo| , 
for a line and parameter range where cc should be at 
its best. In both cases, much better agreement (roughly 
by a factor of 2) is obtained by semiclassical(sc) calcu- 
lations. This means that sc calculations are in fact the 
best available today, in the sense of giving agreement with 
experiment. 

At the foundation of any sophisticated J|,|l6]-|o| sc per- 
turbative calculation is the requirement that unitarity is 
not violated. This is important, of course, since unitar- 
ity violation can lead to a serious overestimation of the 
width 0,^l) . Unitarity is preserved for over 30 years by 
using a criterion, thought to be both necessary and suffi- 
cient. This "fact" has gone unchallenged over this period 
of time. Hence, a minimum impact parameter p m in{v) 
is determined, such that unitarity is satisfied for larger 
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impact parameters, by numerically solving the equation 

(0.1) 

where {. . .} denotes an angular average, V' denotes the 
sc [Q emitter-perturber interaction in the interaction pic- 
ture, a and b denote upper and lower level states respec- 
tively, a' and b 1 denote a complete set of states that per- 
turb a and b respectively, and where d is a number less 
than or equal to 1. The condition d=l is sufficient to 
preserve unitarity, but to keep the expansion parameter 
small, frequently d=0.5 is also used ]l9| . Often, the real 
part of the LHS of the expression in eq.(l) is substituted 
for the absolute value sign, because the imaginary part 
is usually much smaller. The test in eq.l is then carried 
out using the A-function |^,||,|l8| for which analytic ex- 
pressions have been recently given (22). Thus, one hopes 
that perturbation theory is reliable for p > p m i n (v) and 
that this region gives the dominant contribution, so that 
one may either neglect or estimate very roughly the un- 
known (within a sc perturbative framework) p < p m i n (v) 
contribution. In this work, we demonstrate by means of 
a specific case that checking unitarity in this manner 
is not sufficient and the error associated with the break- 
down of this test can be substantial. 

The main problem of the impact approximation has 
always been the so-called strong collisions, i.e. the col- 
lisions at small impact parameters p < p m i ri (v). It has 
been known for a long time that the contribution of such 
collisions can be bounded, but usually some fraction of 
this bound is then used as an additive contribution to the 
width. Whereas for say large scale computations, return- 
ing a single value for the width is desirable, for detailed 
comparisons with experimental results it is best to return 
an error bar for the contribution of these strong collisions, 
which cannot be computed reliably. This is the approach 
adopted here |j| . In other words the theoretical width lies 
between the weak collision width {p > p m in{v)), which is 
presumed to have been reliably computed by the usual 
perturbative (PR) treatment, and the sum of the weak 
and strong {p < p m in{v)) collision widths. Sometimes 
these bounds were satisfactory while other times |7]] 
they were rather too large and reducing them would be 
desirable. 

Although one loosely speaks of small impact parame- 
ters and/or small velocities as giving rise to strong col- 
lisions, it is important to distinguish between collisions 
which fall into this group because unitarity is violated as 
the perturbation expansion fails and collisions which fall 
into this group because the sc approximation is no longer 
valid. Whereas there is no way within the sc approxima- 
tion to avoid the later, the collisions which violate unitar- 
ity are still treatable within the sc impact approximation, 
by means of a non-perturbative(NP) approach. There are 
essentially two such known NP approaches. The first one 
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is an analytic solution under the approximation of ne- 
glecting time-ordering effects, while the second ^3| is a 
fully numerical solution of the Schrodinger equation. To 
avoid ambiguities, we have chosen the second approach 
here. Thus, if the p m in( v ) required to satisfy unitarity 
is significantly larger than then p m in{v) required to pre- 
serve the sc approach, as is very often the case [[| |l0| a 
NP calculation achieves a very significant reduction of 
the error bars given by the PR calculation. 

As a reminder, the half width(HWHM) is written as 

@ 

HWHM = 2nn J pdp J dvv f{v)Q{p,v) (0.2) 

where n is the electron density, f(v) is the Maxwelian 
velocity distribution and Q is given by 

Q(p, «) = {/- S a (p, v)S b \p, v)} (0.3) 

where the subscripts a and b denote the upper and lower 
levels respectively, I is the unit matrix, S is the S-matrix 
and {. . .} denotes angular average. When one formulates 
the problem on the collision axes, one finds for isolated 
lines (jf| 

Q(p, v) = (2J a + l)- 1 J2 (J b m' b lM\J a m' a ) 
(JbmblM\J a m a )[S mb ^ m > b S ma ,m' lL 

-(J a m' a \S c {p, v)\J a m a }(J b m' b \S c (p, v)\J b m b )*} (0.4) 

where the subscript c denotes that the S-matrix has been 
computed for a given direction of the perturber trajectory 
( collision axes). 

For ion lines, the trajectory is parametrized in terms 
of the " time" variable u, defined by f||l8f 

t = s(esinh[u] — u)/v (0-5) 

with s = A7r lornv 2 an< ^ w ith the eccentricity 

All calculations described in this work refer to the sin- 
glet 2s3p-2s3s NeVII 3643. 6A line, which has been the 
subject of a recent experimental study j|. The combi- 
nation of available experimental data and the fact that 
broadening is determined by a small set of three levels 
makes this line suitable for illustrating the important 
points. We first give an example of a pure dipole cal- 
culation. When we calculate Q(p, v) for a velocity of 
2xl0 6 m/sec and an impact parameter of 5A(e — 1.653), 
the PR result is Q(p, v) = 0.05108, and is in good agree- 
ment with the NP result of 0.04997. This is understand- 
able, since Q(p,v) is small(=0. 05108), so the perturba- 
tion expansion should be accurate. But when we attempt 
the same calculation for an impact parameter of 0.68A(e 
=1.0158), which also has a relatively small expansion pa- 
rameter of Q{p,v) =0.16, the NP result is 0.0868. Fig.l, 
which shows the real part of {/ — J7 a (u)J7 b _1 (u)} versus 
the "time" variable u, illustrates what has happened: If 



3 



the U-matrices remain close to the unit matrix, then the 
perturbation expansion is certainely valid. This is the 
case of the solid line, which represents the larger eccen- 
tricity. However, we may get a small Q(p,v) without 
this being the case, as is demonstrated by the dashed 
line, which represents the smaller eccentricity. In other 
words, it is only for small times (or u) that the U-matrices 
evolve according to perturbation theory in the smaller e 
case. Consequently, the unitarity criterion is insufficient 
and the NP result is substantially different from the PR 
one. We note the flat initial and final regions for both 
cases in Fig. I, which show that the U-matrix has indeed 
converged to the S-matrix. 
Fig.l goes here 

Nevertheless, the PR and NP HWHM for a dipole-only 
interaction turn out to be close: 0.535 vs. 0.5094A re- 
spectively. The reason for this agreement can be seen in 
Fig. 2, where the PR (solid line) and NP (dashed line) 
Q(p,v) are shown as a function of p for different veloci- 
ties. Keeping in mind that the average velocity is 2.9xl0 6 
m/sec, it is clear why the differences between the PR and 
NP calculations are small: The main contribution comes 
from weak collisions for which the PR and NP calcula- 
tions are in good agreement. 

Fig. 2 goes here 

Things are different in the quadrupole case where the 
dominant contribution comes from smaller impact pa- 
rameters and velocities. In PR calculations, unitarity re- 
quires a cutoff at larger impact parameters compared to 
a pure dipole calculation, with an associated, often very 
significant, increase in the strong collision error bars. In 
such calculations including quadrupole impact broaden- 
ing often |L8| only the diagonal channels are included. In 
||20| this assumption was discussed and a rough criterion 
was adopted for including the "near hydrogenic" nondi- 
agonal channels. Although tables of the relevant nondi- 
agonal functions exist [|24f| , analytic expressions have not 
been developped for these functions, unlike the dipole 
case, and it would be desirable to do so, in order to be 
able to determine which quadrupole channels should be 
included in the calculation. We here present concrete ex- 
amples where the "near hydrogenic" nondiagonal chan- 
nels must be included: For the line in question, the only 
diagonal quadrupole channel is 3p-3p, since 3s-3s is not 
allowed. However, the 3s-3d channel is allowed and this 
is often neglected |Q. Fig. 3 shows a comparison of pure 
quadrupole calculations. There are basically two regimes, 
the perturbative regime and the strong collision regime 
at small impact parameters. The first thing to note is 
that even in the PR regime, it makes a significant dif- 
ference whether the 3s-3d channel is included or not. On 
the other hand, because perturbation theory is valid, the 
3d-3d channel makes no difference in the results. This 
changes in the strong collision regime and the addition 
of the 3d-3d channel makes an important difference, re- 
sulting in a peak rather than a trough. This result con- 
firms that as we move to smaller p, progressively more 
and more perturbing levels and channels come into play. 
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Fig. 3 goes here 

Fig. 4 is a similar graph, but with the dipole terms 
included. It should be noted that the strong collision 
regime shrinks to smaller p for a pure dipole interaction 
or if the velocity is increased. The remarks made about 
the importance of the "nearly hydrogenic" 3s-3d channels 
still hold true. 

Fig. 4 goes here 

With regard to the effect on the final line profile, the 
electron temperature (average velocity) and density (av- 
erage impact parameter) are very important in determin- 
ing the relative significance of the strong collision regime. 
For the Neon line considered here || , the strong collision 
regime is very important for a dipole plus quadrupole 
calculatiomthe PR result is 0.563Awith an error bar of 
0.477A, as opposed to error bars of about 0.13A for the 
NP results, for which the calculations involving the a) 
3p-3p,b) 3p-3p and 3s-3d, and c) all 3p-3p,3s-3d and 3d- 
3d quadrupole channels yield HWHMs of 0.85, 0.977 and 
0.975A respectively. (The experimental result is ||] 0.865 
A.) It is noteworthy that the 3s-3d channel accounts for 
about 15% of the total (i.e. including the dipole inter- 
action) electronic width, whereas the 3d-3d channel re- 
sults in a decrease of the total electronic width by less 
than 0.5%. These results thus demonstrate the need for 
analytical approximations of the relevant nonhydrogenic 
quadrupole broadening functions p4| , which should then 
be included in existing codes. These results (with their 
error bars) are absolutely rigorous (within the number of 
states used) , if one trusts the "demarkation line" between 
what may and what may not be treated semiclassically, 
and for which the convention in |j| was adopted. It is 
then interesting that even for this case, where relatively 
low partial waves are involved, i.e. under conditions not 
most favorable for sc, good agreement is obtained. Fur- 
thermore, the effects desribed here were responsible for 
the seemingly excess broadening observed in NeVIII ||] 
for which NP calculations give an electron FWHM contri- 
bution between 0.9 and 0.9+0.13 A, with the later num- 
ber corresponding to the (unlikely) maximum possible 
contribution from nonsemiclassical collisions, compared 
with an experimental result of 1.2±0.1 Aat the highest 
density considered. 

In summary, we have demonstrated that the standard 
practice of even the most sophisticated current perturba- 
tive impact theory calculations, cannot enforce unitarity 
reliably and that a more reliable verification of the uni- 
tarity is essential. Such a procedure would involve addi- 
tional information and could be implemented by evalu- 
ating, at intermediate times, the function {/ — UaU^ 1 } 
to ensure that it is monotonic. One possible implemen- 
tation of this could involve using an upper limit of the 
t\ integrations of rather than 00 in Eq.(l); this check 
would allow the diagnosis of the situation shown by the 
dashed line in Fig.l. Even without any additional checks, 
it is still much better to check than not to, since for large 
impact parameters, the unitarity check is reliable. 

Further, it has also been demonstrated that to the ex- 
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tent that quadrupole interactions are important, it is not 
always sufficient to include only the diagonal quadrupole 
channels. An analytic solution to the problem of strong 
collisions is the subject of current research. 
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I. FIGURE CAPTIONS 

FIG.l Real part of {/ - U a U^} vs. u. The solid 
and dashed lines correspond to a NP calculation for v = 
2xl0 6 m/sec and p — 5 and p — 0.68A respectively. 

FIG. 2 Q(p,v) vs. p for different velocities v. The solid 
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and dashed lines represent respectively PR and NP dipole 
calculations. 

FIG. 3 Q(p,v) vs. p for a pure quadrupole interaction. 
&)v= 2xl0 6 m/sec, b)w=3xl0 6 m/sec. The solid, dotted 
and dashed lines correspond respectively to NP calcu- 
lations including all, all except the 3d-3d, and only the 
3p-3p channels respectively. The dash-dotted line cor- 
responds to a PR calculation including only the 3p-3p 
channel. 

FIG. 4 Q(p,v) vs. p for dipole and quadrupole in- 
teractions, a) v=2xl0 6 m/sec, b) v=3xl0 6 m/sec. The 
solid, dotted and dashed lines correspond respectively 
to NP calculations including all dipole channels and, re- 
spectively, all, the 3p-3p and 3s-3d, and the 3p-3p only 
quadrupole channels. The dash-dotted line is a pure 
dipole NP calculation. 
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